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ABSTRACT 
A comparision of the standard and exact methods for taking 
absorption into account has been carried out for the cases 
where Bloch wave degeneracies occur due to the presence of 
non-systematic reflections. Diffraction conditions have been 
chosen so that only three beams need be taken into account 
in the calculations. All the calculations have been carried 
out for materials of differing atomic number. When the 
reflection chosen for dark field image (referred as the 
systematic reflection in this thesis) is in its Bragg 
condition two degenerate Bloch waves are obtained. Under 
these circumstances it has been found that the results 
obtained from the standard and exact methods are in good 
agreement with one another. 

When the systematic reflection is close to its Bragg 
condition, the two Bloch waves are nearly degenerate. Under 
these circumstances marked differences have been found 
between the results obtained from the two methods. These 
differences have been found to increase with increasing 
atomic number. 

In a study of the effect of higher order non-systematic 
reflections it has been found that the differences in the 
results obtained from the two methods increase with 


increasing order of non-systematic reflections. 
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Theselectronsmicrescopemis am. importam.toomunetne 
Stucyr ote (hesstrucitne wot materialise. bineorncers tore the 
microscopist to interpret the experimental micrographs 
obtained from his instrument, it is necessary that reference 
be made to some theory which relates image contrast to 
specimen structure. The interpretation of electron 
micrographs of crystalline materials in terms of electron 
diffraction theories has played an important role in the 
development of electron microscopy as applied to problems in 
material science . Initially, after the development of the 
electron microscope by Knoll] and Ruska: 1932), the 
interpretation of the images obtained from crystalline 
Spee melis Was Gahried GUl im terms) Of Ghe Winemati cal theory 
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kKinematical theory of X-ray diffraction which had been used 
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X-rays and crystals. However a number of authors( for 
examples Blackman) 19395 Von Borries sand Ruska, e19407 
Heidenreich, 1942; Hillier and Baker, 1942: Boersch 
WG42 9194 8e Kinder) 164smsheidenreich andy Ssturkey 9945) 


reported experimental results which could only be partially 
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explained in terms of the Kinematical theory. As a result, 
an interest was developed in the dynamical theory of 
electron diffraction(Bethe 1928). Although this theory will 
be described in detail in chapter 2, for the present it is 
useful to note that Bethe’s approach is quantum mechanical 
in nature and is based on the interaction of fast electrons 
with a crystal as described by the Schrodinger wave 


equation. The wave functions are Bloch waves of the form 


> . F : 1). 
where g is a reciprocal lattice vector, kK is a Bloch wave 


vector and paeic a Bloch wave Fourier coefficient in the the 
direction kK" + ae In solving the Shrodinger equation, the 
lattice potential is expressed as Fourier series which can 


be written as 


where Vg 1S°a Fourier coefficient of ‘theblatticeypotential. 

Bethe’s formulation of the dynamical theory does not 
take the effects of inelastic scattering into account. Some 
of the inelastically scattered electrons will be scattered 
outside the objective aperture and therefore do not 


contribute to the image(see fig.3, chapter 3). These 
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electrons are therefore effectively absorbed and it is in 
this context that the term absorption is used in this 
thesis. As will be shown in chapter 3, absorption can be 
taken into account in the dynamical theory by the 
introduction of an imaginary part of the lattice potential 
(Kambe and Moliére,1970; Dederichs,1972). The lattice 


potential is then complex and can be written as 
V(R) + iV (RB) =D (Vg + iVg )exp(2nig.t) fs3) 
s 


Another important aspect of the theory is the manner in 
which the imaginary Fourier coefficient 1Vg of equation 
(1.3) is incorporated in the theory.The approach widely used 
in the literature is to assume that the imaginary part of 
tne Vat tree* potential ws¥mucn smaller than the real Sant, 
thus permiting first order perturbation theory for non 
degenerate states to be used (Hirsch et. al,1965). This 
method is referred as the standard method in the literature. 
However this form of perturbation theory is only applicable 
to the diffraction conditions for which bloch wave 
degeneracies do not occur (in this thesis the Bloch wave 
degeneracy will be taken to mean an equality in the kinetic 
energies of the Bloch waves (see equation 2.15, chapter 2)). 
In most work performed in the past Bloch wave degeneracies 


were not taken into consideration. As a result the standard 
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method for calculating absorption coefficients, based on 
first order perturbation theory for non-degenerate states, 
has been found to be adequate. During the last several 
years, however there has been considerable interest in 
phenomena, such as the critical voltage effect(Lally et.al 
1972, Andrew and Sheinin 1974A) in which the occurrence of a 
Bloch wave degeneracy is involved. Bloch wave degeneracies 
are also found to occur when non-systematic reflections are 
taken into account in the dynamical theory (Gjonnes and 
Koier 1971; Cann and Sheinin 1972; Sprague and Wilkins 
1970)(for definition of non-systematic reflection see 
section 2.2.1, chapter 2). Sheinin and Cann(1973) and 
Serneels and Gevers(1973) have suggested that perturbation 
theory for doubly degenerate states could be used in 
evaluating absorption coefficients in critical voltage 
calculations. This theory overcomes the difficulty of 
infinite correction to the eigenvectors referred to by 
Spraguesand: Wi lkims’ 41970) y=s Ang@al ternati ves approacheis to 
solve the Schrodinger equation including complex potential 
without making any approximation about the size of the 
perturbation or about how close the kinetic energies of the 
Bloch waves are to one another. This method is referred as 
exact method in this thesis. In this method there are two 


approaches to the solution of the Schrodinger equation. 


These are 
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a. diagonalizing the complex matrix directly (Lally 
Steal 19v 2aethomas #97 2)eand 
b. using generalized perturbation theory(Andrew and 
Sheinin 1975). 
Theesecond approach is more usefull sthan the firs.tesince vit 
gives some insight into the interaction between the Bloch 
waves which results from the introduction of the complex 
lattice potential. 

In a comparison of critical voltage calculations based 
on the standard method with those obtained from the method 
based on exact theory, Lally et.al.(1972) and Andrew and 
Sheinin(1974A) found that errors obtained by using the 
standard method are small! in the lighter elements but ina 
heavy element such as gold, significant errors were 
obtained. Sheinin and Cann(1973) found that significant 
errors could also be obtained in diffracted beam intensity 
calculations in which non-systematic reflections were taken 
into account. A detail investigation of the nature of these 
errors has, however, not been presented in the literature. 
Because the standard theory is so widely used in diffracted 
beam intensity calculations, it was felt important that such 


a study be carried out. The present work was undertaken 
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j1.2 Scope of the present investigation 

In order to carry out an investigation of the errors 
obtained in using the standard theory for the case when 
non-systematic reflections are taken into account, a 
comparision of calculations based on standard method with 
those based on exact method has been carried out: 

deme Ol, Oli Pdeld One Concditd ONS TOW tchmed ib loch wave 
degeneracy occurs in the presence of a 
non-systematic reflection, 

DP ee OCC nimactiony cond lth Ons ton whi che (Won. ocheawaves 
are quasi-degenerate(nearly degenerate or having 
close eigenvalues), 

c. for elements of different crystal structure and of 
differing atomic number and 

d. for both low order and higher order non-systematic 


reflections. 
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CHAPTER 2 


ASPECTS OF THE DYNAMICAL THEORY OF 
ELECTRON SDE PRACT DON 


em eee 


The dynamical theory of electron diffraction, as first 
proposed by Bethe(1928) and further developed by Macgillavry 
(1940) and Heidenreich(1949) starts with the Schrodinger 
equation for an electron in a crystal potential Vir}. Since 
the crystal potential V(r) is periodic in nature, it can be 


expressed in the form of a Fourier series as follows 


Ve enh =a Vg exp(2nig.r) 
3 
=(h2/(2me)) ¥ Ugexp(2nig.t } (oma 
3 


where the summation is over all reciprocal lattice vector g. 


iGroniy Melasticescattening is considered, *thercrysirall 


Belen alerscerea™, Vir) = v*() | and therefore, 
Ug Ug 2} 


In addition if the crystal has a centre of symmetry at the 
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Outside the crystal, where V(r) = 0, the solution of the 


Schrodinger equation is a plane wave of the form 
wiP) = exp(2niXY.r) (2ta) 

where magnitude of the wave vector % is such that 
(h2 *)7 2m =e eE 6285.) 


where E is the potential through which the electron was 
accelerated before entering the crystal. In the crystal the 
electron wave functions are Bloch waves, that is plane waves 
modulated by a function which has the periodicity of the 
latttvees (MAccordingly the solutuems tosthe Schrodinger 
equation are of the form of equation(1.1}(see chapter 1). 

If the expression for V(r) and wi?) are substituted in the 
Schrodinger equation, a set of equations is obtained which 


can be written as 
i meee (1) 
(iii sani daceem leu yal Caner > H0 (2.8) 
h 


where prime on the summation indicates that the term h = 0 


is excluded and 
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where 
a iE (2.7) 
— > (tv) — ¢¢ 
Ka oe (2.8) 


K is the magnitude of the electron wavevector in the crystal 
deem connec lingsforumean icrysta liipotenti alee winnts 
homogenous set of equations gives the general relationship 
between the Bloch wave Fourier coefficients com the Fourier 
coefficients of lattice potential Ug and the Bloch wave 
vectors k‘”. These equations(2.6) are referred as the 
‘dispersion equations’ by Bethe(1928). 

Now the question which arises at this point is what 
approach should be taken in solving the set of equations 
(2.6) . One approach which have been extensively used in the 
past(Whelan and Hirsch 1957; N.Kato 1952 ) is to assume that 
only two beams are important. Under these circumstances the 
set of equations (2.6) consists of only two equations for 
wiicheana lytical solutionsvexists thowever sin thespnesent 
work we are primarily interested in studying the effects of 
non-systematic reflections(see section 2.2.1) and therefore 
more than two beams must be taken into account. If N 
diffracted beams are considered then (2.6) consists of N 
equations for which ana lytical solutions! genendiiyaecanmnat 


be found. Numerical solutions of these equations must 
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therefore be adopted. There have been a number of 
formulations of the multi-beam dynamical theory which can be 
used to obtain such numerical solutions (see for example 
Sturkey(1957), Fujimoto (1959)). The formulation employed 
in the present work is based on expressing (2.6) in the form 
Cian eigenvalue equation(Howie and Whelan 1961). This 
formulation is particularly useful since any number of beams 
can be included and the form of the equation leads itself to 
straightforward numerical solution. 

The manner in which the eigenvalue equation is derived 
can be understood by considering Fig.1 . The values oe 
satisfying (2.6) lie on a surface called the dispersion 
surface. The wave points of the Bloch wave excited inside 
the crystal are obtained from the condition that the 
tangential components of the Bloch wave vectors must be 
equal to that of the incident beam. Thus a line is drawn 
through the point T (in fig.1) and perpendicular to the 
erystal surface to intersect the branches of the dispersion 
surface. In the symmetrical Laue case this line will be 
parallel to the Brillouin zone boundary. Now let y'” be the 
distance between possible wavepoints w'! and the point T and 
let Sg be the distance of reciprocal lattice point g from 


the Ewald sphere(both quantities are measured in the z 


direction i.e., normal to g). It can be seen from fig.1 


that 
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Ly \2 Qs 
Ke © Tors Kau es sy cost, + 1 costa) )2 
= K2 - ( K + Xo Sq) COSBg )2 
Smee Vie EsalicosOad! thet mit tic Sq)/2K} cos8g | 
(2.9) 
If it is assumed that K>>g and k>>(¥- sg} then 
cosGq = 1 eee, Oe 
and 
(ics sy Wack WE= #0 (D411) 
Equation (2.9) then reduces to 
. (0) wig ot (1) _ 
K 2 ran) tse Shy Sg ee ee) 


and Bethe’s equation (2.6) can be written in the eigenvalue 


form 


’ we) 
isp. 
Here C is a column vector whose elements Ca are the 


Fourier coefficients in equation (1.1)(see chapter 1) and A 
jsmcamathix witheelementswA,, = 0; Aggu-aSqeancma qua 
Ug-,/2K where g #h . The eigenvalue Y'’ is related to the z 
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Co ster (2.14) 


If N beams are considered the matrix A will be an NX N 
matrix. Given such an input matrix, computer programs can 
be written to calculate the N eigenvalues and associated 
eigenvectors | om iet@ad ae eet Le .BByacarnyingeoutesuch 
matrix diagonalization procedures for a series of crystal 
orientations, it is possible to determine the shape of the 
different branches of the dispersion surface . In order to. 
be ani to refer to the different Bloch waves excited, it is 
convenient to label the Bloch waves and their corresponding 
branches of dispersion surface. The method generally 
adopted in the literature is to use the integer labels 
Wea tind: sewed has in order of the decreasing values of ve 


(Humphreys and Fisher, 1971 ). 
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excitation coefficients and determine the extent to which a 
particular Bloch wave is excited in the crystal. In order 
to calculate the diffracted beam intensity Ig , equation 
(2.15) is sorted out into components in the direction of Gg. 
If these components are then multiplied by the phase term 
exp(-2niK.r), the intensity of the beam g at depth z in the 
enystais isefoundiitorbe 


(L) (4) 


Paizie- | O, (z) meet Boot Cy exp(2ni Y'’z) | 2 ee en, 


where 0g(z) is the amplitude of the beam g. 
The excitation coeffients in equation (2.16) can be 
calculated from the boundary conditiors at the top surface 


one the -ceysital*atime: 


0,(0) = 1, Q,(0) = 0 ig #0) (2.17) 


Ci eu eS 


(ty, 
where C is a square matrix having the elements Ca in the 
g-th row and i-th column, X is a column vector containing 
the excitation coefficients x in the i-th row and Bp iS el 


column vector containing the diffracted beam amplitudes 
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WrOy0 ve. "ate the tops sumrace of "the: crysta lt mmine excitation 
coefficients in equation (2.16) can then be found by 
obtaining numerical solutions to the non homogenous set of 
Inearwequalions (2518 ihe etfort slomsol Ves thicuse lor 
equations can be greatly reduced when only elastic 
scattering is considered by noting that the A matrix of 
equation (2.13) in this case is real and symmetric. Hence 
the matrix C of normalized eigenvectors is orthogonal and 
ou = C’ , where C’ is the transpose of Cong] ttequa tione2 ane 


is multiplied from the left by © , then 
X= Cou (2219) 


and it can be immediately seen that the excitation 
GOCHPICTENts x" are given by the elements in the first row 


Oem Ten 
oe ete. (2.20) 


The expression for the diffracted beam intensity therefore 


becomes 


pate m= Wyo exmianii(e’ 2 ial? fOR2 1h 
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where ve eben Thus it can be seen that the contribution 
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of a particular Bloch wave to the g-th diffracted beam 
amplitude is determine by the magnitude of b, which is 


reffered as the excitation amplitude in this thesis. 


2.2.1 Systematic and Non-systematic Reflections. 

When a high energy electron beam is incident on a 
crystal it is generally found that a number of different 
sets of planes are close to satisfying their Bragg 
conditions. In forming conventional strong beam images the 
crystal is usually oriented so that only one low order 
diffracted beam is strongly excited. If this reflection has 
a reciprocal lattice vector or then all reflections 
corresponding to the vectors ng where n is an integer are 
called systematic reflections(Hoerni, 1956). Similarly 
those reflections which have reciprocal lattice vector, fh 
not colinear with fo} are called non-systematic reflections. 
Whether a particular systematic reflection is termed 
systematic or non-systematic depends entirely upon the 
low-order diffracted beam being used for imaging purposes. 
Thus a reflection which is termed systematic in one 
eqituation may Upom choosing sa) diifenentslow-ordereren lection 
for study, be termed non-systematic. This is illustrated in 
fig(2) which shows a computed electron dittracti1ontpattern 
fon ae al OrientatlOnminedecrySta lew ithe CcesEnUC tlhe -aumia 
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image formation, then the reflections (004), (006), (002), 
(004) and (006) are designated as systematic and all other 
reflections such as (131), (133) are designated as 


non-systematic. 


The dynamical theory as it was originally formulated 
did not take into account relativistic effects. These 
effects have been considered by Fujiwara (1962), who 
developec a relativistic form of the dynamical theory using 
the Dirac wave equation. His results show that the 
non-relativistic theory can be corrected for relativistic 
effects by employing two simple corrections. These consist 
of replacing the non-relativistic wave length by the 
reliativistically corrected one and multiplying the Fourier 
coefficient terms, Ug by a relativistic mass correction 


term 
B= (1 = v2/e2 } (2.23) 


where v and c are the speeds of the electron and light 
respectivly. 

[hemimagimdnyepdnurchm tne |OUR Teh = coc nici eiiismo maune 
lattice potential Ug (in case of inelastic scattering , see 


chapter 3) is relativistically corrected by multiplying Ug 
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by v7'B (Howie, 1962).The validity of these corrections has 
been confirmed experimentally by Hashimoto(1964), Dupouy et 


ale (1965) andéGoringeret fal.1966): 
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CHAPTER 3 


PPREG@USeOR INELASTIC SCATTERING 


Cee la Lroduction 

The theory described in the previous chapter does not 
take into account the effects of inelastic scattering. 
Electrons which are inelastically scattered outside the 
objective aperture do not contribute to the image and can 
therefore be considered to be absorbed (see fig.3). The 
effects of absorption are taken into account in the 
dynamical theory of electron diffraction by the addition of 
an imaginary part to the lattice potential (Dederichs( 1972) 
or Kambe and Moliére (1970)). There are several approaches 
which have been adopted in the literature to the solution of 
the Schrodinger equation including the imaginary potential 
(see chapter 1). The two approaches (i.e., standard method 
and exact method) which have been used in the present study, 


Wiliebe discussed in this chapter. 


3.2 The Standard Theory for Taking Absorption into Account 
In the standard method it is assumed that V’ (r)<< V(r), 


therefore permitting first order perturbation theory for the 


non-degenerate states to be used. The effect of the 
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A schematic diagram showing the paths of electrons scattered 


inelastically outside the objective aperture (as for example 


AA'). 
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perturbation iV’ (Ff) is to change the energy of the Bloch 


wave by an amount e AE given by 


Gay 


eAE = i(¢b ate ie sac (3.4) 


This energy change can be related to a change Ak, in the z 


component of the Bloch wave vector, i.e, 


b 
»~ 
"1 


2 -(me/h2K) AE = iq C322) 
where qu is the absorption coefficient of Bloch wave bem 
This coefficient can be calculated from the integral in 
equation 3.1 which gives the expression 


= L 
oe ae 


ee EP 
imeem 2 226g ca si | (3. 3a) 
In matrix form this equation becomes 
gilt oxmne oe tng: 6 (3. 3b) 


where A’ is a matrix contaning the elemento U aan) hele ene 

g-th row and h-th column. It can be seen from equation 2.14 
that the change iq’ in er wi! produce =aes imninl ams changes 
aie Accordingly, when the effects of inelastic scattering 


are included in the dynamical theory, the expression for the 
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diffracted beam intensity given in equation 2.16 becomes 
tid eet) i) t 
Pra) le Os(z) Wee LEC, Ca OxXpi2o iy Sic al | eam 


By examining this relation it can be seen that the 
diffracted beam intensity is effectively attenuated or 
absorbed in the crystal due to the presence of the 

exp (-2nq z ) terms . Calculations of the Bloch wave 
absorption parameters indicate that some Bloch waves are 
more highly absorbed than others. This has been explained 
by several authors(for example, Hashimoto, Howie and Whelan, 
1962) in terms of Bloch wave channeling and gives rise to 


the well Known effects of anomalous absorption 


3.2.1 Anomalous Absorption 

The differences between the absorption coefficients, q(t) 
gives rise to effects commonly referred to as anomalous 
eUSOupLIOn el W@eCrCenmLOMmlUS Lhd Ven WeSeRe Gee ic manramis 
convenient to examine the variation of the diffracted beam 
intensity with thickness in the perfect wedge-shaped 

crystal. Consider first a situation where absorption is 
neglected and where it is assumed that there are only two 
Bloch waves with equal excitation amplitudes Cou Fale sa 


case , diffracted beam intensity is found to vary 


sinusodially with crystal thickness as shown in fig(4a). 
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Fig. 4 
The variation of diffracted beam intensity with thickness 
when (a) no absorption (b) normal absorption and (c) 


anomalous absorption is taken into account. 
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Fig.4(b) shows the same situation except that both important 
Bloch waves have been assumed to have equal absorption 
coefficients. It can be seen in fig.4(b) that the mean 
intensity decreases with increasing thickness but the 
sinusodial variations are still present. This is termed as 
normal absorption and is characterized by the fact that the 
visibility or contrast of the fringes as measured by 
Re. Tl inin ad) Aibimos Pel me edoesenotachangeaAwithadeothwain 
the crystal. Fig.4(c), on the other hand, shows the 
intensity profile obtained when one Bloch wave is absorbed 
much more strongly than the other. This profile illustrates 
the effect of anomalous absorption in a wedge-shaped 
crystal. In thicker regions the intensity variations have 
essentially disappeared due to the strong attenuation of one 
of the Bloch waves. However considerable diffracted beam 
intensity is still present due to the weakly absorbed beam. 
Besides having an important effect on intensity profiles in 
perfect crystals, anomalous absorption has also been found 
foulead tomimportant eitects an thewsimages otecrysiall 


defects. 


3.3 Exact method for Taking Absorption into Account 
The exact method, based on the generalized perturbation 
theory, has been presented in detail in a paper by Andrew 


and Sheinin(1975), and therefore in this section this theory 
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is only summarized. 
When the effects of absorption are taken into account 
by introducing an imaginary part to the lattice potential, 


the eigenvalue equation (2.13) becomes 
=a PGC es) 


where iA’ is a matrix in which the elements Ay, are given 
Dyethesreldtion Aqua Uaen /2ke elheamatr i xmAniA sini sano 
longer Hermitian and thus has complex eigenvalues and 
eigenvectors. The real eigenvalues and eigenvectors Y" and 
cM in equation (2.13) have therefore been rewritten in 


10m) (7?) 
anda Ge 


equation (3.5) as the complex quantities Y 
respectively. In the generalized perturbation theory 
approach to the solution of equation (3.5) the perturbed 
eigenvalues Y'™” can be found without approximation by 


solving the eigenvalue equation 


Beast Sev uncer (3.6) 


where B is a matrix which can written as 


meh) : 49) 3 
Y'+i¢q ish Me ad ot 
ey 
ae agze (2), epee iq’ 24 Hew (3.6a) 
; 1) (32 (3) (3) 
felt i qikaeegy( tad Get tem 


owt Gh b@eoP Bia tt oa OR ee Wee Se ae See ee 


ee hal LA 
Tad ae? aT 


Govtd 8° 24 esrancle Ane ron ft ‘iniah @ af. cn. 
on. at ' ‘20s satan al Ro dagle 3 a" rh tae 
os equbsVnegts (oignha 2Eer marl HAs age nse 

unig | ¥ g86'aewhaars ora '2e Buf avira is fkad ef, 2" 
ni aval 2) es", oad eae ayn Hal a) et aaa 

a > Ons VV eer}? 1ARUD staat ans a; i210) oT a > 
gens doniaiwisg- ast feagnsy saint. yfantts 

bert era art $3.4)! bo Seuhs oo nol tulce att oF ABO 1° 
NG Hot Thairso ies /lvariww he? od seo * * yaniv 

hor faves albeit tala an? ¢ 


w) 


(ee) | ; ng OMY a “or 
o8 setaston Neo ral, Sita se © 


“. 


31 


t) 


The diagonal elements Y'"+ iq' On theeareethe sumo 
eigenvalues of the unperturbed equation 2.13 and absorption 
coeffiecients calculated by standard method, while the 


values of off-diagonal elements ape are given by 
ty) th) 6) 
Cumeeee=! 11-21 Gab Mm] Eb atime 
S) 
2 amie Came (307 a) 


and qe” = qt? ete) 


The off-diagonal elements gt) represent the mixing which 
occurs between a particular pair of unperturbed Bloch waves, 


iheseigqenvector cr” can be found from the relation 
Oar Cnakaa , (3.8) 


where C is the matrix in which the columns are the 
eigenvectors of the unperturbed equation (2eee and 2 is 
an eigenvector of equation (3.6). If C is orthogonal and 
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fee veenorma lized a mitienec will also be normalized. 
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The expression for the diffracted beam intensity in the 


exact theory is similar to the one in the standard theory 
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32 
(see equation 3.4) and can be written as 


#(n) 


in) 
Cy SxoOv nie | az ia 2 (3.9) 


The eigenvalues Nedaan this equation can be found by 
diagonalizing the matrix B and Bloch wave Fourier 
coefficients CamecaD be found by substituting the 
eigenvectors or Ayand Bamatrnicessintovequal jon is 18) 8. ihe 
excitation coefficients lean be found by numerical 


solution of the non-homogenous system of linear equations 
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CHAPTER 4 


CeO ene OA ECU eiman Ss 


Hic dauisiere Ve @ vols 

The work presented in this thesis relies extensively on 
numerical calculations based on the standard and exact 
theories presented in the last two chapters. It is 
therefore important to indicate the procedure used to 
Gad leu late andydi agonal ize the indirices A see equaltiongzs ts) 
Hic eeeseescdUaciOheo..ojmand then tOscaleulatesabsor pion 
coeffiecients gM, excitation amplitudes ion and 
diffracted beam intensities. These procedures will be 
discussed in this chapter. Since the work presented in this 
thesis involves Bloch wave degeneracies, it is also 
jmpOOrtame tO indicate the diftraction cond) Gionsea te wWomeh 


the degeneracies are obtained. This will also be discussed 


PiieuiSeChiaOler. 


4.2 Setting up the A matrix 
ieee Ca leu ationeor sune Di adOnalme ements Gue., 


Matrix. 
The diagonal elements of the A matrix (see equation 


2.13, chapter 2) can be expressed in terms of the deviation 
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PUG fs) 
The interesection of the Ewald sphere with the zero order 
Laue zone. The orientation is specified by the tie point 
method. The point L is the perpendicular projection of the 
centre of the Ewald sphere (tie-point) on the zero order 


Laue plane. 


Figure 5 
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eis@ha 18) 
The dispersion surface and the Ewald sphere construction 
Showing parameters required for calculation of diagonal 


elements of the Asmat x. 


37 


Brillouin Zone — >» 


' 

---,-------------------- f----—-qz— Branch i of Dispersion Surface 
{ 
! 


Tie-Point 


Sphere of Radius K 
Centered at O 


Zero Order 
Laue Zone 


Figure 6 


e5abhu® doimaqgad tc 1 ecesa 


iniaTeiT — 


We; 


M awiheS igh wige 
f) io Perne. 


‘ { 


unt) get 
wrens ved \ / j 


' | | ; 
- Pade 2. | 
a” 2 ———— | is be eine? aang hlawS | 


Et a me mw. ee -s=—Aam& > 2 ne 


e motel 


: 


2405 ows veh? iin J 


eewgs 


38 


parameter s, where S), is defined to be the distance from the 
reciprocal lattice point h to the Ewald sphere in the 
z-direction (see fig.6). When the position of the Ewald 
sphere is specified using the tie-point method (in this 
method the crystal orientation is specified by the 
co-ordinates of the point L which is the perpendicular 
projection of the tie point in to the zero order Laue plane 
bseent i gr Seanceii0.,o Gee then values of svacay be calculated 
with the aid of Fig(6). In this figure, the symmetrical 
Laue case has been assumed and the point L is the progee tien 
ofethe £tue-pointMigonithe zero’ ordertiaues pianeser 1t can) be 
seen that the value of Sy for some reflection h in the n th 


order Laue zone is 


Using a Taylor series expansion of the two square roots and 


retaining terms up to second order, gives 


Senor OK Scrat sare ae Kala, = ace ae (4.2) 


This method of calculating the diagonal elements of the A 


matrix will be referred to as the s-value method. 
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It is of interest to note that if only the first order 
term of equation 4.2 is retained, the s-values for 


reflection g in the zero order zone becomes, 


Spee olh/ 2K lise Grave tu) (4.3a) 
=plat/ 2K Myths om) Wt GeH 
SA ig Cee =o 5) 
= (1/Kigiz,- 9/2 ) 
= gle (4. 3b) 


where A#® is defined in fig.(6). This first order 


approximation is used by many authors(see for example Head 


Eitenaslee, bhSias 


The A Matrix. 
The off-diagonal elements of the A matrix(see section 
Do weform a) crystal icomposedaof .only sonestype ofeatom are 


given by 
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fF(sin6,_, /X) is the atomic scattering factor for electrons 
and B is the Debye Waller factor for the material 
considered.” The scattering factors used in this thesis were 
calculated by Doyle and Turner(1968) using relativistic 
Hartree-Fock atomic wave fuctions. Values for the 
Debye-Waller B factor are also given in literature(Ibers and 


Vainshtein Coes 


4nseSetting up the B matrix 

Ine@aragomalmelementssot the Bematri xd scesequationmos.c) 
are obtaimed by first diagonalizing the real matrix A\see 
section 4.4). The absorption coefficients of the standard 
theory (see section 4.8) and the off-diagonal elements gq 


can then be calculated from the equation (3.7a). 


inorder 3vo obtain Y'° and Csi the matrix A(see 
equation 2.13), which is real and symmetric, was 
diagonalized by using the subroutines EHOUSS, EQRT2S, EHOBKS 
in IMSL subroutine pakage(1979). It should be noted that 
there are other subroutines in the IMSL(1979) which could 
haverbeen wusedem.) combarisconeo; al) thesroutinesswas 
carried out and it was found that the subroutines used in 
the present study require the least computing time and the 


results are more accurate than those obtained by using other 
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subroutines(Cann 1973) 

To calculate Y and C the complex matrix Paice 
equation 3.6) was diagonalized by using the IMSL( 1979) 
SUDnoOuUtIi ne =e LGCG. 
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In this thesis three beam calculations have been 
Carried out with the three reflections (000), (220), (133)" 
In a comparison of experimental results with many beam and 
three beam dynamical theory, Cann and Sheinin(1974) found 
that there are three important Bloch waves when (133) 
reflection is close to its Bragg condition. In order to 
verify that no other reflections have a significant effect 
on the final result, both three and multi-beam calculations 
of Bloch wave parameters were carried out for different 
materials which have been considered in this thesis. Good 
agreement was obtained between the results of both sets of 
Calculations, andicating hatwein=thtcoecasemuhess moles tirnee 
beam approximation is adequate for predicting the effects of 
a Bloch wave degeneracy in the presence of non-systematic 
reflections. In order to determine the effect of higher 
order non-systematic reflections a study has been caried out 
with the higher order non-systematic reflection (pce 

One question which arises when performing the 


multi-beam calculations is how many reflections should be 
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A computed diffraction pattern for a [334] orientation in a 
crystal with fcc structure showing the 11 reflections which 


were considered in the many beam calculations. 
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included. In view of the fact that the computing time 
required for a particular calculation increases considerably 
with the number of reflections taken into account, it was 
necessary to develop a criterion which selects only those 
beams which are relevent to the results obtained. The 
procedure used in the present work was to increase the 
number of reflections included until a stage was reached 
when the addition of further reflections did not change the 
intensity profile or the values of Bloch wave parameters 

Ga Da 


reflections were considered for the multi-beam calculations 


and qin any significant way. In this way eleven 


(see fig.7). 


4.6 Diffraction conditions for which the Bloch wave 


degeneracies are obtained 


The conditions for the occurance of a Bloch wave 
degeneracy in the presence of a non-systematic reflection h 


(Cann and Sheinin 1974; Gjonnes and Hoier 1971) are given by 
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where S4 and s, are the deviation parameters of reflection g 
and h (see section 4.2.1) and U, and U, are parameters 
related to the Fourier coefficients of the lattice 
potential(see equation 2.1, chapter 2). For the particular 


case of (220) and (133) the expressions (4.5) and (4.6) 


reduce to 
Sc = tae Sa) 
7 (yee wi 
a pe Ett) Va oaul 
2K Ua20 
ancucamiarlyerCmmethe casesot 91000) smo 20) mui ods) 
Seay =o) (4.5b) 
Uja7 — Ure. 
ee ee TS Care (4. 6b) 
137 ZK U220 


In these particular cases two Bloch waves will be degenerate 
only when (220) reflection is in its Bragg condition and 
deviation parameter of (133) and (137) are given by 

Edud pions 4 sOduchGe4nODPmmn i Sa1S i VUSIha ecm Gece wiiien 
shows that Bloch waves 2 and 3 are degenerate at AGj3; = 
-0.056;,, ( 48;3, is the deviation of (133) reflection from 


ioe etaod Cond 1 tion wedi Vehemst dC GlOnS#O? sc-asmaic (133) 


Bragg angle).This value of A8;3; as taken from the curves 
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PUG (3 
The variation of ove en eee) eEWaLe | A873 as given by a 
three beam calculation including (000), (220). and (133) 


reflections. AB 559 = 0.0 


Rae hey we) 
The: Variation. of yt = ibsesn) iraliels A8733 as given by a 
three beam calculation including the (000), (220), and 


(1338) "reflections. Ac = 7-0. 16 
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corresponds to a deviation parameter S33 =-0,107¥%100° A in 
agreement with the value obtained from equation (4.6a). 

When the (220) reflection is not in the Bragg condition a 
degeneracy does not occur. However when A6j20is small, 

Bloch waves 2 and 3 are nearly degenerate, a condition which 
will be referred as quasi-degeneracy. This is illustrated 

in fig.9, which shows the quasi-degeneracy at AS; = -0.07674; 
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ae ecalculationsoi Absorption .coefticientsssexci tation 
amplitude and Diffracted beam intensity using the Standard 
me thod 

The absorption coefficients q' were calculated by using 
equation (3.3a)(see chapter 3). In order to carry out this 
calculation the imaginary component of the lattice potential 
Ug must first be obtained(see equation 4.7). This was done 
by using the ratios Ug /Ug given by Humphreys and Hirsch 
(1968). The excitation amplitude oe | as given by 
ct cl” (see equation 2.21), can be obtained from the 
eigenvectors of matrix A (see section 4.4). The diffracted 


beam intensity Ig was obtained from equation 3.4 (see 


chapter 3) 
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4.8 Calculation of Absorption coefficients, Excitation 
amplitude and Diffracted beam Intensity using the exact 
me thod 

The absorption coefficients gq’? in exact method were 
obtained by diagonalizing the complex matrix B where the 
imaginary part of the eigenvalue gives the aegenet ier 
coefficient(see chapter 3,section 3.3). The excitation 


(t) Se) 
amp 1itudes 1%, |, given by XC.” | where the Bloch wave 


q 
Fourier coefficients ot found by diagonalizing the matrix 
B(see section (4.4) and the excitation coefficients Naeescnd 
by solving the linear equation (2.18)(in this case C ,X and 
UMarcecomp 16x eUSTnOMLnem IMSL (9/3 Meubroutimes BE Q2 CRs estae 
diffracted beam intensities were obtained by writing a 
program based on equation (3.4) and in which ae aH A 


and gq’ were calculated by using exact method. 
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CHARTER 5 


Re OUME SP ANDEDISCUSSION 


DRaeli eROGUC ton 

The Bloch wave absorption coefficients gs and 
excitation amplitudes toe | have been calculated as a 
function of AB;y33 ( S8;j33 is the deviation of the reflection 
(133) from its Bragg condition) by using the standard and 
exact methods. Calculations have been carried out for both 
three beam (including (000), (220), and (133) reflections) 
and eleven beam cases(see fig.7). The intensity profiles 
(j.e., diffracted beam intensity plotted as function of 
crystal depth) at the point where Bloch waves 2 and 3 are 
degenerate ( s,.,= 0 } or quasi-degenerate (| s,,,| is close 
to zero} have also been obtained by using both methods. 
Calculations have been carried out for three materials 
macludingsSigeCuvand AUPand ancompardsongot laiinthe 
calculations based on the standard method with those based 


) 


on exact method has been made. In the plots of gq‘ and 


(W) , 
| ay | as a function of A8+3, , Bloch wave one has not been 
220 133 


included. The principle reason for this is that the values 


(1) (1) 


Grace dnc| Vee | are not effected by degeneracy and 


therefore both the standard and the exact methods give the 
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same results for this Bloch wave. A study has also been 
carried out to determine the effects of a higher order 
non-systematic reflection (137). Results of all these 
studies will be presented in the subsquent sections of this 


chapter. 


o.2 A Comparison of Absorption coefficients, Excitation 


Cee a ea eine ae Se eee eee 


Standard and Exact methods. 


5.2.1 The case of degenerate diffraction conditions. 


It was indicated in the previous chapter(see section 
4.6) that in the presence of the non-systematic reflection 
(133), a Bloch wave degeneracy occurs when the (220) 
FeTleCction 1Ssine1ts. Bragg conditions. ihe first 
calculations to be presented are for the medium atomic 


number material copper at an accelerating voltage 100Kv. 
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Figures 10(a,b) and i1t(a,b) show the variation of | o.., | 


) 


and qi’ with AGyz33 obtained by using standard and exact 


methods respectively . From these figures it can be seen 


that there is no significant difference between the values 


(UES) 5 


Of o | or g°? obtained from either method. A point of 


220 

interest which should be noted, however is that there is an 
Qu) " 

interchange in the values of | Call and Can for Bloch waves 

2 and 3 which occurs at the point of degeneracy. The value 


of 48.,, at the point where the interchange occurs has 
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Pichia ONCE 0) 
The variation of excitation amplitude be (i302 5s )ewath 


TONE Clee a can = 0 calculated by (a) the standard 


Searle 220 
and (b) the exact methods. 
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Fils Wey) 
The variation of Bloch wave absorption coefficients Aare 
with \8733 FO Ge Clliemd v A8noq = 0, calculated by using (a) 


the standard and (b) the exact methods. 
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LAB Lee 


Percentage differences in the absorption coefficients and 
diffracted beam intensities calculated by using both 
standard and exact methods at the point where Bloch waves 2 
and 3 are degenerate. Three beam calculations were used and 


the results shown are for silicon, copper and gold. 


Nore, = 0 Accelerating voltage = 100 kv. 
Crystal thickness for Cu and Au = 1500 A 


Crystal thickness for Si = 15,000 A 
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different values for the two methods. However, because this 
interchange occurs at a particular value of 46;3;, the 
values of the parameters in equation (3.4)(see chapter 3) 
for i =2,3 remain unaffected on passing through the 
degeneracy point. As a result there should be no 
significant difference in the intensity profiles obtained by 
the two methods in this case, as indicated by the data in 
table(1). 

Similar calculations have also been carried out for 
gold and silicon. The results obtained with these two 
elements were found to be similar in character to those 


obtained for copper (see table 1) 


5.2.2 Results in the case of quasi-degenerate 
Cillania CullOomcOnGul Ons 

In section 4.6(chapter 4) it was indicated that when 
there is a small deviation of (220) reflection from its 
Beagousecond itlonn Gb 1CCh Waves s2ndiGecmcalsbemsCOnS IOCheCmtOnDe 
quasi-degenerate (see fig. 9). In this thesis a study of 
the quasi-degenerate state has been carried out for a value 
Ope AGen., = 0.00018. > 5. 

The results obtained for copper are shown in 
pic caro manger Ce icv CUGeS= 12 lar lo eeshOWstae 
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(i = 2,3) with 48;., obtained 
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Pika, WAL ey) 
The variation of (a) excitation amplitude loss? Ge 268) 
and (b) the variation of aoe Cp eae) wale A987 33: 
Calculations were carried out for Cu using both exact and 


standard methods, at a value of A®# = -0.00019 
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Palle) VES 
The variation of diffracted beam intensity with crystal 
thickness for Cu obtained by using the standard and exact 
methods. Calculations were carried out at the point where 


Bloch waves 2 and 3 are quasi-degenerate, i.e. at A8733 = 


-0.050- 


133 and A@ = -0.00018 : 
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these Figures that in the case of the standard method the 
interchange in the values of g® ele || ed occur over a 
range of values of A833 | i.e., the interchange is 
continuous). In the case of the exact method the 
interchange is similar in character to that obtained in the 
case of degenerate diffraction conditions(not shown in 
fiche @)eeeerOMmnr tom) aeb ee Gecansbesseccomthavusa teinespoint 
where Bloch waves 2 and 3 are quasi-degenerate, i.e., at 


A65,, = ~0.053336;23, the standard method gives 


(2> €3) 


229 | = | $,5,|1> whereas the exact method 


values are quite different. As a result it would be 
expected that the intensity profiles obtained with the two 
methods would be quite different as shown in figure (13). 
It should be noted from this figure that there is no 
appreciable change in the overall shape of the intensity 
profiles obtained from the two theories and the peak to peak 
Spacing 1s MOt molicably altered. “ihew two; theories® 
however, give significant differences in diffracted beam 
intensity. 

Similar calculations have been carried out for silicon 
and gold. The results obtained were found to be similar in 
character to those for Cu as illustrated in figures (14), 


(a Seand stable n\ 2). 
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Pita. al4 
The variation of diffracted beam intensity with crystal 
thickness for Si obtained by using the standard and exact 
methods. Calculations were carried out at the point where 
Bloch waves 2 and 3 are quasi-degenerate, i.e. at A873 = 


and Aé@ ==), 00013 : 


Sue 220 220 


E33 


uns IRs) 
The variation of diffracted beam intensity with crystal 
thickness for Au obtained by using the standard and exact 
methods. Calculations were carried out at the point where 
Bloch waves 2 and 3 are quasi-degenerate, i.e. at A@= 


iocien 


-0.115653. and AG 509 = -0.0001655,- 
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Figure 14 
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TABLE 2 


Percentage differences in absorption coefficients and 
diffracted beam intensities calculated by using standard and 
exact method at the point where Bloch waves 2 and 3 are 
quasi-degenerate. 3 beam and 11 beam calculations were used 


and the results shown are for Si, Cu and Au. 


Mere = SO sOes ely. Accelerating voltage = 100kv 
Crystal thickness for Cu and Au = 1500 A 
Crystals thickness for oi. = 15.000 By 
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FAG Ode) 
GO) 


The variation of the Bloch wave absorption coefficients q 


with A85 33 for Cu at AS 509 


using both standard and exact methods. 11 beams were taken 


= -0.00016 559; calculated by 


into account tn) the ‘calculations. 


bag. 16(b) 
The variation of diffracted beam intensity with crystal 
thickness for Cu obtained by using the standard and exact 
methods. Calculations were carried out at the point where 
Bloch waves 2 and 3 are quasi-degenerate, i.e. at A873 = 
-0 03673, and A859 = -0.0001655,- 11 beams were taken into 


account inetiesca bCull ati oOnse 
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E — Exact 
S - Standard 
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Multi-beam calculations have also been performed in 
order to find out if the three beam approximation is 
adequate. Calculations have been carried out for the three 
elements (i.e, Cu, Au and Si) for quasi-degenerate 
diffraction’ conditions (ie, when AGlo5 =9- 0 0001Gss0) 

Eleven reflections have been considered (see section 4.7) 
for these calculations. It has been found that the results 
obtained are similar in character to those obtained in three 
beam case, as illustrated by the results obtained for Cu are 
shown in fig.16(a,b) and figures (12,13) (see page 48,49) 
and table(2). 


De2eoeet fects of higher orden snon2zsystematic 
reflection. 

A study has been carried out to see if the effects 
found in the previous sections are affected by the order of 
the non-systematic reflection. The reflection considered 
for this study is (137) and calculations have been carried 
out for both Si and Cu. Both degenerate ( A6230= 0 ) and 
quasi-degenerate ( A8jz25= -0.00018220) diffraction 
conditions have been considered. The first calculations 
have been carried out for Cu. It has been found that in 
case of degenerate diffraction conditions, the results are 


similar in character to those obtained for (733) reflection 
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Fiepe Wey) 
The variation of the Bloch wave absorption coefficients 
ae with A873 for Cu. Calculations were carried out by 
using both standard and exact methods at a value of 


Aé =~-0 00018 . 
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Panes UTE bles) 
The variation of diffracted beam intensity with crystal 
thickness for Cu obtained by using the standard and exact 
methods. Calculations were carried out at the point where 
Bloch waves 2 and 3 are quasi-degenerate, i.e. at A= 


aon 


-0.03673, and AB 509 = -0.00016,5)- 
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Figure 17 (a) 
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TABLES 


Percentage differences in absorption coefficient and 
diffracted beam intensities calculated by using standard and 
exact method at the point where Bloch waves 2 and 3 are 
degenerate and quasi-degenerate. 3 beam calculations were 
used including (000), (220) and (137) reflections and the 


results shown are for Si and Cu. 


Crystal thickness for Cu = 1500 A 
Crystal thickness for Si = 15,000 A 


Accelerating voltage = 100kv 
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and therefore no significant differences were found in the 
results obtained with standard and exact methods (see 
table(3)). In the case of quasi-degenerate diffraction 
conditions, however, significant differences have been found 
in the results obtained with the two methods as shown in 
figures 17(a,b) and table(3).Similar calculations have been 
carried out for Si. The results obtained are similar in 
character to those obtained for Cu (see table 3). 

It is of interest to compare the results obtained for 
the low order reflection (133) with those of higher order 
reflection (137). In order to carry out this comparison 
percentage differences in diffracted beam intensity have 
been calculated for different crystal thickness as a 
constant value of A6,,, (see table 4). Calculations of 
diffracted beam intensity using the two theories have also 
been carried out for different values of AB za9 as a 
constant crystal thickness (see table 5). It can be seen 
from the table(4) and table(5) that the differences in the 
results obtained by using both methods are greater in the 
case of (137) reflection than in the case of (133) 


reflection. The reason for this will be discussed in the 


Next —Ssection: 


shi nl onic? ew -dnsol Vingte tn 
wee) oboriien ebrate “itw & 

nat icant in etommagebetanin 19 sens ant. 
waa naed. svat ears Hi Tneat Fi note: pomey: 
ai mente ge eburiten sted dtrw t . 
fend aver! andtteldeles “efimié>(tistds? one ta, 
a¢ aif imhe ote mats lo giivecs of! «(ld eh 
(C wiowt. ase? iy 4g? banthyas Sea ets 

ae? Gsntetde alivesr ant silos of ‘deweatni io et 
asecd réioi Yo begil Ayiw Sil) polgasttes seh wat “9 | 
naetesgens 2frid tuo VAG ed oaoce PLETELY notiowht 
Wet el bereiat tise bimwivy ) o? aaunewetith egs insted: 
«ss ae aaenaatett \6fagit> Ins wi! thse? Getatuol soe 

to asctiafuel 4G ih e@PRSE RaD), GegB A te sui sv INE? 
cain avad serve? owt ath oe ae VGiesatot mad tele 
5S get % esvisv-thetett(> wh ive bet 

rage sa rao 11 «4 )}¢ etd gen) esenmaln? leleves or 
ea? mii Asomeisttit ony tant leleldet ome (kietees aft mot 
att nt Sieg x6 chodtemortod offen yd tented | 
(S65) Fo seen ark mf want ap hhamt tes wet io 

sit nf Baseunel® sd tw Bhd? 36% a, inet toeh te" 


tan 


73 


TABLE 4 


Percentage differences in diffracted beam intensities 
calculated by using standard and exact method at the point 
where Bloch waves 2 and 3 are quasi-degenerate. 3 beam 


calculations were used and the results shown are for Cu. 


AB555 = ~0.00016225 Accelerating voltage = 100kv 
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TABLE 5 


Percentage differences in diffracted beam intensities 
calculated by using standard and exact method at the point 
where Bloch waves 2 and 3 are quasi-degenerate. 3 beam 


calculations were used and the results shown are for Cu. 


Accelerating voltage = 100kv and Crystal thickness for 
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ome DASCUSSTION 

The results described in section 5.2.1 showed that for 
degenerate diffraction conditions i.e., when the systematic 
reflection is in its Bragg condition, there are no 
significant differences in the absorption coefficients and 
diffracted beam intensities obtained from the standard and 
exact methods. The results described in sections 5.2.2 and 
5.2.4 on the other. hand showed that for quasi-degenerate 
diffraction conditions i.e., when the systematic reflection 
is clese to its Bragg condition there are significant 
differences in the results obtained from the two methods. 
The manner in which these differences arise differ in three 
important respects. The first of these can be seen from the 
numerical calculations in table(5) (see page 59) which show 
that the differences occur over a range of values of 

Re en il) Seetablemshowse ina te themoignini canta dit crences 
in the diffracted beam intensities obtained when the two 
methods for calculating the absorption coefficient were 
used: voccun! Over@a range tofavalues of AC LSP irompabout 
“00000 Geesitols Ota Oseoninicaseton) ecierer lection: 

The second important feature of the results presented 
in this thesis can be seen by comparing the results for 
different materials as shown in table(2}) (see page 52). 
These results show that the differences in absorption 


coefficients and diffracted beam intensities obtained from 
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the two methods are greater in the case of Au than in the 
case of Cu, which indicates that the errors increase with 
increasing atomic number. From tables (1), (2) and (3) (see 
bages™ 4698 52eand'56) Git@candbevseenaithat@the results 
ebtaimed tor diamond cubic crystal (Si) Parelsimivanein 
character to those obtained for the face centered cubic 
6GyStals  GCueanceAune 

The third important feature of the results presented in 
this thesis can be seen by comparing the results obtained 
for low and higher order non-systematic reflections in 
tabies (4) and (5) (see page 58, 59). These results show 
that the differences obtained from using the two methods 
increase significantly with order of the non-systematic 
reflection. The underlying reasons for this behaviour wil 


be discussed in the next section. 
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the systematic reflection is close to its Braqg condition. 
The results described in section 5.2.4 and tables (4) 

and (5) (see pages 58,59) showed that the differences in the 

results obtained from the two methods for calculating 

absorption coefficients are greater in case of (137) 

reflection than for (133) reflection. This can be explained 


in terms of the characteristics of the B matrix. The most 
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TABLE 6 


B matrices for three beam calculations in which (a) (000), 
(220), (133) and (b) (000), (220), (137) reflections are 


inc luded 
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important point to note is that the amount by which the 
standard theory eigenvalues (M4 iqule which are equal to 
the diagonal elements of the B matrix, are modified by the 
off-diagonal elements depends on the value of these elements 
as well as the differences between the diagonal elements 
(Andrew and Sheinin 1975). It can be seen from table(6) 
that the differences between the diagonal elements of the B 
matrix are less for (137) reflection than in the case of 
(133) reflection. It is clear therefore that the difference 
in the absorption coefficients obtained from standard and 
exact (diagonalizing the B matrix) methods should be greater 
inweasemor  iicmimret lec tionmethansina the casetotmticc) 
reflection as has in fact been found (see the right column 
Ofe tabdess). 

The reason for the decrease in difference between the 
eigenvalues ea. with increasing order of non-systematic 
reflection was explained by Cann(1973). In a determination 
of extinction distance in the presence of higher order 
non-systematic reflections Cann found that the minimum 
difference between the eigenvalues is equal to U,/V2k. This 
shows that the difference is directly proportional to 
Fourierecoefficient of the lattice potential ofethe 
non-systematic reflection involved. Since these 
coefficients in general decrease in magnitude with 


increasing order of non-systematic reflection, the decrease 
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in differences between the eigenvalues with increasing order 


of the non-systematic reflection is explained. 


5.4 Conclusions 

The original purpose of the work undertaken in this 
thesis was to determine the errors obtained from using the 
standard theory to calculate absorption coefficient in the 
dynamical theory of electron diffraction for the case where 
Bloch wave degeneracies occur due to the presence of 
non-systematic reflections. In order to carry out such 
investigation a comparison of the calculations based on 
standard theory with those based on exact theory has been 
Ganmied Outi. 

In attempting to assess the significance of the results 
in this thesis the first point to consider is the importance 
of taking non-systematic reflections into account in 
dvitpactron icontras tecaicuhat honsmarinminesodas Hemost 
calculations of diffraction contrast have been based on 
either the two-beam dynamical theory or the multi-beam 
theory taking systematic reflections into account. In 
practice, the assumption that only systematic reflections 
are present is open to question since non-systematic 
reflections must to some extent be excited. In the 
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reflections the electron microscopist is confronted with 
problem of choosing either standard method or the exact 
method to calculate absorption coefficients. 

In this thesis it has been shown that there are no 
significant errors obtained in using the standard theory to 
calculate absorption coefficient when diffraction conditions 
are such that a Bloch wave degeneracy is obtained in the 
presence of non-systematic reflection, significant errors, 
however, are obtained in the case of quasi-degenerate 
diffraction conditions. It is important to appreciate that 
in practice it may be very difficult for the electron 
microscopist to judge, from the diffraction pattern, wether 
degenerate or quasi-degenerate diffraction conditions 
ebtaine —In this situation significant “errors, mayebe 
obtained by using standard method to calculate absorption 
coefficient and therefore it is important to use the exact 
method. It has also been indicated in this thesis that the 
errors obtained by using the standard method increase with 
increasing atomic number and with increasing order of the 
OMe sVStemabic nef leccliona sl fll S tnenerOrescleGanpeatnd tatnere 
are many possible situations for which the exact theory 
should be used when non-systematic reflections are taken 
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single non-systematic reflection is taken into account. In 
practice there may be a number of non-systematic reflections 
excited which may have significant effect on the results 
obtained. Since each non-systematic reflection can give 
rise to a degeneracy, it is possible that the errors 
described in this thesis may be considerably greater in 
practice. Further investigation is of course required to 
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